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We construct a relativistic potential quark model of D, Ds, B, and Bs mesons in which the 
light quark motion is described by the Dirac equation with a scalar-vector interaction and the 
heavy quark is considered a local source of the gluon field. The effective interquark interaction is 

described by a combination of the perturbative one-gluon exchange potential Vcoui(r) = —(,/r and 
the long-range Lorentz-scalar and Lorentz-vector linear potentials S'i.r.(»') = (1 — X){ar + Vo) and 
Vi.i.{r) = X{ar + VQ), where ^ A < 1/2. Within the quasiclassical approximation, we obtain simple 
asymptotic formulas for the energy and mass spectra and for the mean radii of D, Dg, B, and Bs 
mesons, which ensure a high accuracy of calculations even for states with the radial quantum number 
rir ~ 1. We show that the fine structure of P-wave states in heavy-light mesons is primarily sensitive 
to the choice of two parameters: the strong-coupling constant as and the coefficient A of mixing 
of the long-range scalar and vector potentials 5i.r.(r') and Vi.r.(r). The quasiclassical formulas for 
asymptotic coefficients of wave function at zero and infinity are obtained. 
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I. INTRODUCTION 

The heavy-light quark-antiquark (Qq) systems, be- 
ing QCD analogues of relativistic hydrogen-like atoms, 
are ideal objects for investigations and permit very pre- 
cise experimental verification of quantum theory results. 
Theoretical description of the mass spectra and decay 
probabilities of such composite objects requires the con- 
struction of a consistent theory of bound states, which 
should be based on the fundamental principles of local 
quantiim field theory and use its apparatus [1]. However, 
calculating these characteristics of composite systems di- 
rectly in the local quantum field theory is not always pos- 
sible, because the only known calculation method in this 
theory is still based on the perturbation theory, while the 
nature of creation of a bound state of interacting parti- 
cles must undoubtedly be determined by nonperturbative 
effects. 

When constructing the theory of bound states, the 
most effective way of leaving the framework of the per- 
turbation theory is the use of dynamic equations. The 
point is that even if wc c;an construct kernels of dynami- 
cal equations only in the lower orders of the perturbation 
theory, developing methods for solving them exactly or 
approximately (but without using th(^ pc;rturbation the- 
ory) allows taking into account nonperturbative effects of 
interaction when evaluating observable characteristics of 
the bound states. In a nonrclativistic case, such a theory 
is formulated in the language of the classical potential 
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using the dynamical Schrodinger equation. But at large 
bound energies, the corresponding theory must be essen- 
tially relativistic. In this regard, the way to solve this 
problem was indicated about half a century ago based 
on using the dynamical equations in the local quantum 
field theory, examples of which are the Bethe-Salpether 
equation [2], the quasipotential equation [3], and other 
equations [4]. 

The Dirac equation with a mixed scalar-vector inter- 
action plays an important role in the contemporary de- 
velopment of the relativistic theory of bound states. It is 
valuable because it provides an adequate mathematical 
model for a wide circle of problems in hadronic physics 
in which it is possible to transit consistently from a 
two-particle problem to the external field approximation. 
This equation indicates the presence of the spin and spin 
moment for the quark and antiquark, and the problems 
of describing fine and superfine structures in the energy 
spectra of heavy-light {Qq) mesons, which arc the QCD 
analogues of hydrogen-like atoms, arise naturally from 
this equation. Treating the Dirac equation in the limit 
of an infinitely heavy quark Q as an equation for a single 
light antiquark q (similarly for the case of hydrogen-like 
atoms), we can study several important aspects of the 
theory of heavy-light quark-antiquark systems, in partic- 
ular, the relativistic dynamics of the light antiquark q 
in the external field of the heavy quark Q, the Lorentz- 
structure of the long-range component of the Qq inter- 
action, the fine structure of the spectrum of heavy-light 
mesons, and the influence of the spontaneous breaking of 
chiral symmetry on the spectrum, etc. 

The mathematical theory of the Dirac equation with 
a scalar-vector interaction was developed in [5] (see [6- 
8] for a detailed bibliography). Certain progress was 
achieved in constructing the exact solutions of equations 
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of this type with potentials corresponding to different 
types of interaction [5] . However, in most cases attempts 
to construct exact solutions of this equation for more or 
less realistic potentials encounter difficulties that have 
not yet been overcome. The known methods of investi- 
gation of this equation approximately (the perturbation 
theory in the coupling constant, etc.) do not provide 
complete knowledge about the behavior of the wave func- 
tions and mass spectrum in the most interesting domain 
of values of the coupling constant for hadronic systems 
containing one light quark together with one {D and B 
mesons) or two (doubly heavy S and fl baryons) heavy 
(anti) quarks; relativistic and nonperturbative effects evi- 
dently play an important role in such systems. Therefore, 
when constructing approximate methods for investigat- 
ing bound states of the Dirac equation, nonperturbative 
methods, in which the expansion parameter in the po- 
tential is not considered small, are especially important. 
One of the most widely used methods is the method of 
asymptotic expansion in the Planck constant h, which is 
called the quasiclassical approximation. 

The apparatus of quasiclassical asymptotic behavior 
which has arisen soon after making a quantum mechan- 
ics is now one of the most powerful and (frequently) also 
a unique way of study of a wide class of problems of 
theoretical and mathematical physics. Fundamental na- 
ture of a quasiclassical method is caused, in particular, 
by its embodiment of the important (for physical theo- 
ries) principle of correspondence of quantum and classical 
problems. 

In the one-dimensional spectral problems of a quan- 
tum mechanics the strict mathematical formulation of 
this principle composes the Wentzel-Kramers-Brillouin 
approach (WKB). Historically first of all, namely the one- 
dimensional WKB approximation began to be applied for 
solving some of the quantum mechanical problems and 
in problems of the collisions theory. The WKB method 
has allowed to obtain general expressions for reflectivi- 
ties and transition of particles through a potential bar- 
rier, for quasiclassical phases and scattering amplitudes 
in a central field, for Bohr-Sommerfeld quantization con- 
ditions, for probabilities of nonadiabatic transitions at 
pseudocrossing of levels. 

The rigorous theory of quasiclassical asymptotic ex- 
pansions including the scattering problem together with 
spectral problems, was constructed in Maslov's funda- 
mental monograph [9] and subsequent articles [10]. The 
WKB method for fermions satisfying the Dirac equation 
with a purely vector interaction (including states lying 
near the boundary of the lower continuum) was devel- 
oped in detail in [11-13]. No wonder that the framework 
of the quasiclassical methods has solely allowed to obtain 
many approved results in the known theory of superheavy 
atoms [14]. 

The construction of quasiclassical solutions of the 
spinor equation with a scalar-vector interaction was re- 
ported in [15, 16]. The scheme of quasiclassical quantiza- 
tion proposed in [16] allows to make clear the connection 
of quasiclassical asymptotic behavior in spectral prob- 



lems for the Dirac equation in external scalar and vector 
fields with the Lorentz structure of interaction potentials 
corresponding to them. 

To elucidate how the Lorentz structure of interaction 
potentials is manifested at solving spectral problems of 
hadron physics, in the present article we study the be- 
havior of a relativistic spin- 1/2 particle in the presence 
of both the scalar and the vector external fields with 
potentials of confining type by WKB method. More def- 
initely, at creating of the relativistic version of the poten- 
tial model taking into account the Lorentz structure of 
potentials of interquark interaction, we relied on Cornell 
model offered in [17] in which the effective color Coulomb 
attraction at small distances r and string interaction at 
large r are considered. Corresponding to such model of 
Q^-interaction the effective potential (EP) U{r,E) of the 
Dirac equation complicatedly depends on energy E and 
on mixing coefficient A (0 ^ A ^ 1) of scalar and vector 
long-range potentials and has essentially different form 
at A < 1/2 [18], A > 1/2 [19] and A = 1/2. For all fur- 
ther consideration, the important moment here is that 
the same model with a scalar-vector variant of interac- 
tions leads (depending on range of values of the mixing 
parameter A) both to bound states and to decaying (qua- 
sistationary) ones. Therefore, other field of application 
of the used potential model concerns to rather impor- 
tant problem related to study of quasistationary states 
of quantum objects. For this reason as an example of 
the chosen model of interaction, in the present article 
the problems of theory of decaying states are considered 
together with spectral problems. 

The structure of the article is as follows. In the next 
section, we set out the construction of the quasiclassical 
solutions of the Dirac equation with a scalar-vector cou- 
pling. In Sec. Ill, we consider the model of the interac- 
tion of a relativistic particle with both scalar and vector 
fields given by the potentials V{r) ~ —£./r + \{ar + Vc))) 
S{r) = (1 — A) (or -|- Vq) and study the dependence of 
the effective potential on the Lorentz structure of the ex- 
ternal field. In Sees. IV and V, in the framework of the 
considered model with mixing parameter ^ A < 1/2 
we obtain simple asymptotic formulas for the energy and 
mass spectra and for the mean radii of D, Ds, B, and Bg 
mesons. Our results are compared with results of numer- 
ical calculations and existing experimental data. Sec. VI 
is devoted to deriving asymptotic coefficients of the wave 
function at the zero and infinity in quasiclassical approx- 
imation. In Sec. VII we calculate the energy spectrum of 
the massless fermion in the external scalar field described 
by the combined funnel potential. 



II. QUASICLASSICAL APPROXIMATION FOR 
THE DIRAC EQUATION WITH A VECTOR AND 
SCALAR INTERACTION POTENTIAL 

After separation of angular variables in the Dirac equa- 
tion with a mixed scalar- vector coupling (c = 1), system 
of equation for the radial wave functions F{r) and G{r) 
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is of the form 



dF k 

h— + -F- [{E - Vir)) + (m + 5(r))] G = 0, 
ar r 

h^--G+ [{E - V{r)) - (m + S{r))] F = 0. 
dr r 



(1) 



Hereafter, we use the notation F{r) = rf{r) and G{r) — 
rg{r), where /(r) and g{r) are the radial functions for 
the respective upper and lower components of the Dirac 
bispinor [20] , E and m are the total energy and rest mass 
of the particle, S{r) is the Lorentz-scalar potential, and 
the potential V{r) up to a multiplier coincides with the 
zeroth (temporal) component of the four- vector potential 
Afj, = (Ao, A), where A = 0, V{r) = -eAo{r), and e > 0. 
In system (1), k — hk, where the quantum number k = 
tO + 1/2) for ? = j =F 1/2, j is the total angular moment 
of the fermion, and I is the orbital moment (for the upper 
component of F{r)), and hence \k\ = j + 1/2 = 1, 2, . . .. 

In [16] the system (1) was consecutively solved using 
the known technique of left and right eigenvectors of the 
homogeneous system. For the effective potential (EP) of 
the barrier type (see Fig. 1) 



F — l/^ 

U{r,E)^-V + S+^—^ 
m Zm 



(2) 



quasiclassical expressions were obtained for the wave 
functions in three regions. 

i) in the classically allowed region tq < r < ri the 
asymptotic approximation of the WKB radial functions 
F and G are oscillatory: 



Fir) - C± 
G(r) = G^sgnk 
where 



E -V + m + S 



1/2 



p{r) 

E-V-m-S 
p(r) 



cos 9i, 

1/2 



(3) 



cos 82 



p{r) = [{E - F(r))2 - (m + S{r)f - (k/r) 



U{r,E) 



21 1/2 



(4) 




ri 



FIG. 1. The form of the EP U{r,E) of the barrier type; ro, ri, 
and r'2 are roots of the equation = 0. 



is the quasiclassical momentum for the radial motion of 
a particle, and 
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1 / V + S' 1 

2 \m + S - E + V ^ r 



the normalization constants Cj related to the 
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states with fc > and fc < are denoted by the su- 
perscripts + and — . The normalization constant is 
determined by the relation 



\ct\ = 



E-V{r) 
p{r) 



-1/2 



dr 



1/2 



(5) 



where the quantity T coincides with the period of radial 
oscillations of a classical relativistic particle with the en- 
ergy E in the potential well < r < ri . 

ii) in the below-barrier region ri < r < r2, the 
q uantity p takes purely imaginary val ues, p — iq, q — 
yj{m + S{r)Y - (e - V{r)y + {k/r^. Oscillating-type 
WKB solution (3) is then continued here by a solution 
exponentially decreasing as the distance increases in the 
classically forbidden region ri < r < r2. For states with 
fc > 0, we have 



cxp 



{m + S)V' + {E-V)S' 



2qQ 



dr 



-Q 



m+S -E+V 

(6) 

where Q = q + \k\r^^. The solution for the states with 
fc < can be found in [16]. 

iii) in the "exterior" classically allowed region r > r2, 
a quasistationary state is associated with the divergent 



VpP 



exp 



ip 



im + S)V' + {E ~-V)S' 



2pP 



dr 



iP 



m+S -E+V 



(7) 



This expression must be used to study the states with 
fc > (fc < see [16]); here, P = p + i\k\r~^ and the 
radial momentum p{r) is again positive. 

The quasiclassical representations (3)-(7) constructed 
are invalid in small neighbourhoods of turning points rj 
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(j = 0, 1, 2). To bypass these points and match the sohi- 
tions one can use the Zwaan method [21] that allows to 
establish relations between normalization constants: 



3 ^2 



± r 



X exp ■ 



E-V{ri)^m + S{ri) 
(771 + 5) V + {E-V) S'' 

2gg _ 



dr 



(8) 



Also in [16] the following quantization condition, de- 
termining the energy (position) of the boimd state E in 
the mixture of the scalar and vector potentials, was ob- 
tained: 



quasiclassical calculation methods developed in present 
article will be useful also for problems of hadronic spec- 
troscopy in the range of small quantum numbers. 

The spectral problem for the Dirac equation with the 
potentials S{r) and V{r) of the confining type that con- 
sidered in subsequent sections illustrates applying these 
methods to problems in hadronic physics. Other types 
of the potentials S{r) and V{r) and also a more detailed 
mathematical description of the WKB method for the 
Dirac equation with a scalar-vector interaction can be 
found in [16]. 



III. THE DEPENDENCE OF THE EP U{r,E) 
ON THE LORENTZ STRUCTURE OF THE 
EXTERNAL FIELD 



kw 
pr 



dr 



(9) 



Here, 77,. = 0, 1, 2, . . . is the radial quantum number. The 
new quantization rule (9) differs from the standard Bohr- 
Sommerfeld quantization condition [22] by the relativistic 
expression for the momentum p{r) and by the correction 
proportional to w{r), which takes into account the spin- 
orbital interaction and results in the splitting of levels 
with different signs of the quantum number k. 

Having calculated the flux of the particles outgoing to 
infinity by means of quasiclassical formulas (7), (8), we 
find the following expression for the level width: 



r= ^exp[-20]. 



(10) 



r = 2 



^,r. n = /(,-^5;|*. (U) 

b 



which is valid for both signs of k [16]. 

The obtained quasiclassical formula (10) is the rela- 
tivistic generalization of the well-known Gamow formula 
for the width of a quasistationary level. The nontrivial 
moment of such a generalisation is the modification of 
expression for the period of oscillations T and the oc- 
currence of the additional factor in the preexponent of 
expression (10) that depends on a sign of the quantum 
number k and is caused by the spin-orbit coupling in the 
mixture of the scalar S{r) and vector V{r) potentials. 

The described scheme of quasiclassical quantization 
is applicable for deriving the asymptotic behaviors of 
level energy and corresponding eigenfunctions, generally 
speaking, at "large" quantum numbers. However, in well- 
known exactly solvable spectroscopic problems of the rel- 
ativistic quantum mechanics, for example for scalar S{r) 
and vector V{r) potentials of Coulomb {S{r) = —^'/r, 



Vir) 



-i/r) or oscillatory {S{r) = V{r) — 
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w > 0) types, formulas of quasiclassical quantization re- 
produce an energy spectrum precisely (even for lowest 
states) [IG]. This circumstance allows to suppose that 



The simplest model of the interaction of a relativis- 
tic spin-1/2 particle simultaneously with both scalar and 
vector external fields, which we meet below when calcu- 
lating the quasiclassical spectrum of relativistic bound 
states (see Sec. IV), is governed by the potentials 

Vir) = ycoui(r) + ViAr) = - ^ + Mr), (12) 
Sir) = Si.r.ir) = (1 - X)vir), ''vir) = ar + Vo, 

where Vq is a real constant, ^ is the Coulomb coefficient, 
and A is the parameter of mixing between the vector and 
scalar long-range potentials Vi.r.(r) and S'i.r.(r); < A < 
1. Below in this section, we do not restrict the value or 
even the sign of the parameter a. 

The relation between the EP ?7(r, E) and initial poten- 
tials (12) directly entering the Dirac equation is rather 
complicated: Uir,E) depends not only on r and model 
parameters (12) but also on the level energy E and on 
the total moment j. What is especially important for 
us here is that the EP C/(r, E) takes essentially different 
forms for the cases A < 1/2, A > 1/2, and A = 1/2. 

Our goal is to investigate the behavior of the EP 
Uir,E) at large and small r. Substituting Vir) and 5(r) 
of form (12) in (2) and keeping only the most singular 
terms when r — ^ only the leading terms (in r) when 
r — )• 00, we obtain: 

(1 - 2X)a^ 



Uir,E) 



E 



2m 
f 777 



■r' + . 



X^ 



1 



2777 

7^ 



ar + . . 



00, 



0, 7^ = fc2 _ 



(13a) 
(13b) 

2 ^ . (13c) 

- 2777 

For rather large values of the Coulomb coupling constant 
5 in EP C/(r, E) at small distances r the singular attrac- 
tion oc arises, which can lead to "falling to center" 
known in quantum mechanics [22-24]. In order to show 
this, we consider the finite trial function which is nonzero 
in the range < r < tq. In accordance with the Heisen- 
berg's indeterminacy relation (j?^)ro > 1/4, whence fol- 
lows 



2m 



2mr^ 



ro < 1. 
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At ^ > |fc| = j + 1/2 the spectrum of the effective Hamil- 
tonian H is unbounded below, because (H) — )■ — oo at 
vq — > 0. In the classical mechanics such situation cor- 
responds to a particle falling to a force centre. In the 
case of the Dirac equation the energy eigenvalues En^k 
have the squared- root singularity at ^ — >■ |fc| that leads 
to states with complex energy E^^k at the formal con- 
tinuation of the one-particle solutions into the domain 
£, > \k\ — j + 1/2. It does not mean, however, that 
at ^ > the Dirac equation has no solutions. For 
determination of energy levels in this case it is neces- 
sary to impose some boundary condition at r — >■ (that 
is equivalent to determination of the self-conjugate ex- 
pansion of the energy operator [25]). From the physical 
point of view the formulation of a boundary condition at 
zero means the cut-off of the (color-) Coulomb potential 
VcouK?") at small distances, i.e. taking into account finite 
sizes of an extended source (heavy quark, nucleus): 



V{r) = 



— Xv(r) for r>rN, 
r 

-—/[ — ) for < r < Tat 



(14) 



(here tn is the cut-off radius, and /(O) < oo). The form 
of the cut-off function f{r/r]^) is determined by distribu- 
tion of (color) electrical charge over volume of a nucleus 
(quark) (see [26, 27]). 

Under the requirement ^ |cr|~^/^ which is satisfied 
for small values of the Coulomb parameter ^ < level 
energy En^k depends weakly on a concrete form of V{r) 
at small r. Therefore, we carry out the passage to the 
limit (rjv — >■ 0) of the pointlike Coulomb potential in 
(14). If ^ > |fc| then the dependence of En^k on the form 
of cut-off function f{x) becomes rather strong which is 
characteristic for all problems with "falling to center" . 

Let us now consider behavior of EP U{r,E) and wave 
functions at large distances r in more detail. First note 
that only the quadratic term {S^ — V'^)/2m is essential in 
the asymptotic domain in formula (2) for A 7^ 1/2 and has 
the behavior (1 — 2A)(T^r^/2m when r — > 00. It is hence 
obvious that for any sign of the parameter a, the EP 
U{r, E) of model (12) under consideration (at sufficiently 
large distances) is an attractive potential for A > 1/2 and 
a repulsive potential for A < 1/2. Both types of behavior 
(i.e., attraction for A > 1/2 and repulsion for A < 1/2) 
are purely relativistic effects related to the fact that the 
interaction of the fermion with the scalar external field 
S{r) is added to the scalar quantity m, the particle mass, 
while the vector potential V{r) is introduced into the free 
Dirac equation minimally as the temporal component of 
the Lorentz-vector A^. 

It is clear from what was said above that for A < 1/2, 
the EP [/(r, E) of model (12) is an unboundedly increas- 
ing (as r increases) confining potential with only a dis- 
crete spectrum of energy levels; it is then essential that 
the quadratic dependence of the EP U{r,E) on r (and 
hence the confinement property) appears because of the 
relativistic terms (5^ — V'^)/2m. An example form of the 
EP U{r,E) for A < 1/2 is shown in Fig. 2. It is amazing 



that bound states are present in composite field (12) un- 
der consideration for A < 1/2 even in the case where the 
initial long-range potential v{r) = ar + Vq corresponds 
to attraction (cr < 0, Vq < 0). 

But for A > 1/2 and an arbitrary value of cr ^ 0, the 
effective Hamiltonian H of the squared Dirac equation in 
external field (12) has complex eigenvalues of energy be- 
cause the EP U{r,E) becomes negative in this case (at 
sufficiently large distances) and less than the effective 
particle energy E = (i?^ — m?)/2ra, which corresponds 
to attraction. Therefore, for A > 1/2, the EP U{r,E) 
of model (12) has the form of a well separated from the 
external domain by a wide potential barrier (for |cr| ^ 1; 
see Fig. 1). It is obvious that the leading contribution 
to forming the barrier of the EP U (r, E) comes from the 
Lorentz-vector component V|.i..(r) of the long-range po- 
tential v{r). Furthermore, as follows from (2) and (15a), 
in the presence of only a vector field (A — 1), the EP 
U{r,E) does not have the confining property even when 
the initial long-range potential v{r) — ar + Vq corre- 
sponds either to attraction (cr < 0, Vq < 0) or to repul- 
sion (ct > 0, Vq > 0). This is the principal difference 
between relativistic potential model (12) under consider- 
ation and the analogous nonrelativistic model in which 
the EP U^ff{r) = -$/r + ar + VQ + l{l + l)/2r^ in the 
radial Schrodinger equation has the barrier for negative 
values of the parameters a and Vg, which results in qua- 
sistationary states with complex energies appearing in- 
stead of discrete levels. On the contrary, if cr > 0, then 
the EP ('') becomes an unboundedly increasing con- 
fining potential with only the discrete spectrum of energy 
levels. The absence of bound states in the Dirac equation 
with a linearly increasing vector potential V{r) was first 
noted in [28] . 

The quasiclassical formulas for the wave functions in 
the domain r > r2 [16] imply the asymptotic form of 
the radial functions F{r) and G(r) as r —> 00. It then 
happens that the wave functions decrease exponentially 
at large distances for ^ A < 1/2 and oscillate if 
1/2 < A ^ 1. As an illustration, we present this asymp- 




FIG. 2. The EP U{r, E) of Dirac system (1) with potential 
(12) in the case where A < 1/2, a > 0, and E > fh; a, b, c 
and d are the quasimomentum roots in (20). 
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totic behavior for the radial function corresponding to 
the upper component of the Dirac bispinor (r — )■ oo): 



exp 



.^lElMr^^ forO<A<l/2, 



exp I I 



\a\ 



(15) 



for 1/2 < A < 1. 



It hence follows that the relativistic solutions for po- 
tential model (12) (depending on the value of the mix- 
ing parameter A) constitute stationary or quasistationary 
systems satisfying different boundary conditions (15) for 
A < 1/2 or A > 1/2. 

We point out one more important particular case re- 
alized at A = 1/2. Substituting potentials (12) with the 
value A = 1/2 in expression (2), we see that the quadratic 
dependence of the "tail" of U (r, E) on r disappears and 
the long-range components V[.r.(^) a-^d 5*1 i..(r) of the first 
two terms dominate EP (2) at large r, which results in a 
practically linear dependence of U (r, E) on r (see (15b)). 
We note that we again obtain a linear confining potential, 
which has only the discrete spectrum, at positive values 
of a, while for negative (sufficiently small) values of a, 
the EP U{r, E) of model (12) has a wide barrier. Because 
of this, level decay by percolation through the potential 
barrier becomes possible, i.e., the bound level becomes 
a quasistationary exponentially decaying state with the 
complex energy E = Er — iV /2. From the analyticity 
standpoint, the above behavior of the EP U (r, E) for 
(T < and (T > allows studying how the discrete spec- 
trum continues from the real axis to the complex plane. 

Summarizing, we can say that varying one of the pa- 
rameters of interaction model (12), the coefficient A of 
mixing the scalar and vector long-range potentials Si,,-, (r) 
and V\,T:,{r), in the interval < A < 1, we obtain qualita- 
tively different forms of the EP t/(r, E): from the confin- 
ing potential with only the discrete spectrum for A < 1/2 
to the potential with the potential barrier and quasista- 
tionary energy levels for A < 1/2 through the physically 
important intermediate case A = 1/2, where the asymp- 
totic behavior (as r — >■ oo) of the "tail" of the EP U (r, E) 
switches from quadratic (15a) to linear (15b) (sec above). 

For the squared Dirac equation in composite field (12), 
the form of the EP becomes more complicated: expres- 
sion (2) for U{r. E) acquires small corrections due to the 
particle spin and the related spin-orbital interaction. It is 
clear from the nature of the conclusions about the behav- 
ior of the EP U{r, E) for A < 1/2, A > 1/2, and A = 1/2 
that the indicated changes of the form of C/(r, E) do not 
change the results qualitatively. 

Everything said above remains valid for the spherically 
symmetric potentials S{r) and V{r) with the powerlike 
or logarithmic behavior (u(r) ^ ar^ , 13 > 0, or v(r) ~ 
^logr) of the long-range part v{r) at infinity. 

Having clarified the qualitative aspects, we now con- 
centrate on a practical application of the above appara- 
tus of quasiclassical asymptotic behavior to heavy-light 
mesons. 



IV. QUASICLASSICAL DESCRIPTION OF THE 
ENERGY SPECTRUM OF HEAVY-LIGHT 
QUARK-ANTIQUARK SYSTEMS 

To use the potential approach to describe properties 
of heavy-light mesons, we must construct the quark- 
antiquark interaction potential. As is known from 
QCD, because of the asymptotic freedom property, the 
Coulomb-type potential of the one-gluon exchange gives 
the leading contribution at small distances (r < 0.25Fm). 

As the distance increases, the long-range confining in- 
teraction (the confinement), whose actual form has not 
yet been established in the QCD framework, prevails. 
The confining potential may have a complicated Lorentz 
structure. For example, it was shown in [29, 30] that the 
interaction of the quark-antiquark pair with a fluctuating 
gluon vacuum field at a finite correlation length results 
in a linearly increasing potential. The spin-dependent 
potential obtained with that approach has a structure 
that is characteristic of scalar confinement. On the other 
hand, the infrared asymptotic behavior of the gluon prop- 
agator of the form D{k^) ~ l/(k^)^ was obtained in [31] 
by analyzing the system of the Schwinger-Dyson equa- 
tions. Such an asymptotic behavior in the static limit 
results in a linearly growing vector confining potential. 
It is therefore most plausible that the confining potential 
comprises a mixture of vector and scalar parts. More- 
over, lattice calculations [32] based on the first princi- 
ples of QCD support a linear confinement proportional 
to r/47ra'(0) (where a'(0) is the slope of the hadronic 
Regge trajectory). From the above considerations, we 
assume that the Qq interaction is a combination of the 
following potentials: 

a. the one-gluon exchange potential Fcoui('') = —i/f, 
where ^ = 4/3as, ag is the strong coupling constant 

as{Q) = 127r/[(33 - 27V/) log(gVA')], (16) 

Nf is the number of quark fiavors, and A = 360 MeV is 

the QCD parameter, 

b. the long-range linear scalar confining potential 
'S'coni ('') =(1 — ^)v{r), where v{r) is determined by ex- 
pression (12), and 

c. the long-range linear vector potential V^onf('') = 
Xv{r). 

The total effective qiiark-antiquark interaction is then 
described by a combination of the perturbative one- 
gluon exchange potential Vcoui(^) and the scalar and vec- 
tor long-range confining potentials S'conf('') and V^onf('')- 
Therefore, the potentials S and V are given by (12), 
where a = 0.18 GeV^ is the string tension, Vq is the con- 
stant of the additive shift of the bond energy, and th(^ co- 
efficient A of mixing between the vector and scalar confin- 
ing potentials is the adjustable parameter, ^ A < 1/2. 
We can consider that the value of a., is approximately the 
same for each family of heavy-light mesons and doubly 
heavy barions and changes in accordance with (16) only 
when we pass from one family to another. 

In the nonrelativistic limit [E m and S,\V\ <C m) 
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EP becomes 



Here, we use the notation 



Ueff « S{r) + V{r) + 



where the sum 



Un.r. = S{r) + V{r) = - ^ + ar + Vb 
r 



(17) 



(18) 



plays the role of an interaction potential 

Therefore, the description of Qq-interaction by means 
of the representation (12) for potentials V{t) and S{'r) 
is nothing else than a generalization of the interquark 
interaction potential (18) to the relativistic case of QCD- 
motivated Cornell model [17]. 

As is known, the potential (18) possesses solely a dis- 
crete spectrum at cr > and is the spherical model of 
the Stark effect in the hydrogen atom at cr < 0. In 
the nonrelativistic problem it is indifferent what of terms 
{V{r) or S{r)) in (18) ensures the confinement of quarks. 
However, the analysis of the Dirac system (1) with po- 
tentials (12) shows (see Sec. HI) that the radial wave 
functions F{r) and G{r) exponentially decreasing (at 
r — >■ oo) and normalised properly can be gained only 
when Sconfir) > Vconfir) at < r < oo. 

We cannot solve Dirac system (1) with potentials (12) 
exactly; we hence use the quasiclassical approximation 
method, which provides a high accuracy even for low- 
lying quantum numbers in the case of scalar and vector 
fields of the Coulomb and oscillatory types [16]. 

Choosing the mixing coefficient in the range ^ A < 
1/2 corresponds to the scalar confinement prevailing. In 
this case, the EP U{r, E) of our model has the form of 
a standard oscillator well with a single minimum (at the 
point rmin ~ and no maximums (see Fig. 2). 

The equation p'^ = 2m.{E — U{r,E)) =0 determining the 
turning points then results in the complete fourth-degree 
algebraic equation r'^ + f + g r"^ + hr + 1 = with the 
coefficients 



2[rh{l- X) + EX] 
(1 - 2A) 

2EC 



£2 



m 



h 



(1-2A) a2 



9 = 



I = 



2^aX 



(1 - 2A) 



(19) 



(1-2A) tT2 



where E = E — AVb, fh = m+ {1 — A)Vo are the charac- 
teristic parameters with the respective meanings of the 
"shifted" energy and the "shifted" mass. This equation 
has four real roots d < c < h < a determined by the 
equalities 



a=-{ + ^(S + A+),6 = -{ + ^(S-A+), 
c = -{-^(H-A-),d = -{-i(S + A_). 



(20) 



f_2g u(2 
4 2, 2,\Z 



1/3 



1/3- 



1/2 



D 



A± = ^jF±—, Z = w + \/-4u3+t;2, 



2 3 3 Iz 



1/3 



1/3 



_ 1 /Z^ 
" 3 U, 

D = -/3 + 4/.g - 8h, u = - 'Sfh + 121, 
v = 2g^ - 9fgh + 27h^ + 27 pi - 12gl. 

For the potentials under consideration, the quasiclassi- 
cal momentum is determined by equalities (4) and (12). 
Using formulas (20), we represent it in the form conve- 
nient for what follows (cr > and ct < 0) 



R{r) 



We integrate in quantization condition (9) over the clas- 
sically allowed domain between the two positive turning 
points ro = b < ri = a, while the other two turning 
points (d < c < 0) are in the nonphysical domain r < 0. 
Using formula (21), we transform quantization integrals 
(9) into the sum of the integrals 



a 

Ji = J p{r)dr 

b 



R 



dr, 



p(r): 



-dr 



dr 



+ 



{r-X+)R J (r-A_)i? 
Lb 6 



dr 



2\a\y/l-2X 
where we introduce the notation 

E + mT\J{E + fhY - 4cr^(l - 2A) 



(22) 



A± = - 



2a (1 - 2A) 



Writing condition (9) in terms of Ji and J2 is advan- 
tageous compared with the initial representation because 
the integrals contained in Ji and J2 can be expressed in 
terms of complete elliptic integrals. 

The particle energy spectrum is determined by quanti- 
zation condition (9), which, after quantization integrals 
(22) are evaluated (see Appendix), becomes the transcen- 
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dental equation 

2y/T^^ (\a\{b-cy 
^/{a-c){b-d) \ ^ 

+iV3n(.,x) + iv,n(^.,x)]+2(i_2,)|,| 

X [(6 - c) (TVsn , x) + N^U {v_ , x)) + N,F (x)]} 



[N^F{x) + N2E{x) 
k 



rir + - I TT, 



(23) 



where F(x), E{x), and n(j/, x) are the complete elliptic 
integrals of the respective first, second, and third kind 
(see formulas (A.l)). The mathematical details of calcu- 
lating integrals of type (22) can be found in [33, 34], and 
the expressions for Xj 3?, and Ni [i = 1,...,7) are 
collected in Appendix because they are rather cumber- 
some. 

Finding an "exact" solution of Eq. (23) in the gen- 
eral case is, of course, impossible, but the situation is 
simplified with the increase in the energy E or in the 
approximation of "weak" long-range field (as compared 
with the Coulomb field). The first case corresponds to 
the fact that for not too large (i.e., for "intermediate") 
values of the parameters and (namely, for ct < 0.2 GeV^ 
and 0.3 < ^ < 0.8), the condition E'^ ^ ct7 is well satis- 
fied for all possible values En^k of the heavy-light meson 
energy levels, and the second case is realized when the 
condition a <C ^fn? is satisfied. In the framework of our 
consideration (i.e., for the physics of heavy-light mesons), 
only the first case is interesting, while the second case is 
most often encountered in approximate calculations of 
those properties of low-lying hadronic states that do not 
depend directly on the presence or absence of confine- 
ment. 

A simple and often effective method for deriving 
asymptotic expansions of integrals of form (22) is to ex- 
pand a quasimomentum p{'r) in a small parameter, the 
interaction, and integrate the obtained series term by 
term. We then indicate two special features of this pro- 
cedure for calculating the integrals Ji and J2 containing 
the small parameter. First, it is obvious from analyzing 
expressions (20) that in addition to the level E = m, 
we must introduce one more characteristic energy level 
E = ffi, which divides the domains of applicability of the 
asymptotic approximations for the quantization integrals 
Ji and J2 obtained below. Using the relations E > m 
and E < fh, we can show that in these two domains of 
the spectrum, the motion is quasiclassical if the condi- 
tion a'y/E'^ <ti 1 is satisfied ioi E < fh and the condition 
aj/E"^ -C 1 is satisfied for E > m. This gives the possi- 
bility of obtaining expressions for Ji and J2 in elementary 
functions using the formal expansion of the quasimomen- 
tum in a power series in a small dimensionless parameter 
(which is aj/E'^ <C 1 or a/^fh'^ < 1). Second, the fur- 
ther analysis depends essentially on the mutual positions 
of the turning points a, b, c, and d. Then, depending 
on the relative values of E and the level fh, we consider 



several typical situations. 

Case A: Let a > and the conditions a <C £,fh^ and 
£^ < m be satisfied. This situation describes deep levels 
whose energy is close to the bottom of the scalar- vector 
well U{r, E). Estimating expressions (20) for the turning 
points in the approximation a j^ffi} ^ 1 and preserving 
only the two first terms in the small parameter expansion, 
we can easily obtain 



E(^-e 

fh-E 



1 - 



E£, 



Vi + 



fh£,T]2 



a 



, Ei-e 



(24) 



fh + E ^ (, 



o- fh-E' (7(1 -2A) fh + E 

Hereafter, we use the notation 



e = ^{Eky-{fh^Y, iJL = Vfh^ - E-^, 
jji = (1 - A)m + XE, r]2 = Xfh + {1- X)E. 



(25) 



It follows from asymptotic expressions (24) that the 
positive turning points a and b depend weakly on a and 
are determined only by the Coulomb field. The other two 
(negative) turning points c and d depend mainly on the 
linear part v{r) of interaction (12), but their values are 
"corrected" by the quantities ^^/(fh^ E), which are due 
to the Coulomb long-range interaction. It is also obvious 
from (24) that for small positive values of a, the turning 
points c and d are sufficiently far from the two points a 
and b and tend to —00 in the limit as ct — >■ 0. 

The properties of deeply lying levels for massive quarks 
(m^^ ct) are mainly determined by the Coulomb po- 
tential. Treating the long-range potential v{r) as a small 
perturbation, we can expand the quasiclassical momen- 
tum p{r) in the domain of the potential well 6 < r < a in 
a series in increasing powers of r/|c| ^ 1 and r/\d\ ^ 1. 
Calculating the table integrals in (9) , whose sum gives the 
value of the quantization integrals Ji and J2 up to terms 
of the order 0((ct/^to^)^), we then obtain the equation, 
which can be easily solved for the level energies. 



Erirk = Eo + XVo + 

(2efhEo 
I Mi " 



2^m2 



V20 



Mo 



O 



k mo 



,(26) 



where Eq = m 



-1/2 



is the Dirac level 



of the energy of the fermion (with the effective mass 
fh = m + {1 — A)Vo) in the Coulomb field, n'^ = rir + 
(1 +sgnA:)/2, and the quantities jj-o, r?io, and 7720 are ob- 
tained from n, rji, and 772 by substituting E for Eq. The 
previously accepted condition ct > is unnecessary here 
because this result remains applicable also in the case of 
negative values of the parameter ct. 

Formula (26) can also be found using the standard per- 
turbation theory, but this involves rather cumbersome 
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calculations. Using quasiclassical formulas (9) and (22) 
dramatically simplifies calculations. As is shown by com- 
paring with the result obtained by numerically integrat- 
ing Eq. (1), formula (26) ensures a good accuracy for 
calculating the spectra of bound systems of heavy quarks 
(for example, QQ mesons; see [35]). 

Calculations, which we omit here, demonstrate that 
in the case A < 1/2 and for (sufficiently small) negative 
values of a, the EP U(r,E) has the shape of a double 
well. If we neglect the barrier penetrability in the region 
c < r < b between the two wells, then the quasiclassi- 
cal quantization conditions in this well can be written 
merely as the conditions on the phase integrals over the 
domain of the quasiclassical motion in each of the wells. 
Quantizing in the left well by formula (9) then results in 
formula (26) above. 

Case B: In the domain E > fh and ct > 0, which is of 
actual importance for the physics of heavy-light mesons, 
a small dimensionless parameter a^/E"^ appears in the 
spectral problem. Imposing the condition a^/E'^ <^ 1, 
we can easily obtain the approximate expressions for the 
turning points from exact formulas (20): 



E-fh 



-E^-e 



E 



E^ 



E + m 



(27) 



E^ 



+ 



ct(1-2A) e 



m 



As can be seen from these formulas, the turning points a 
and b are rather distant from each other, and the above 
expansion for the quasimomentum p{r) is not applica- 
ble in the whole integration domain. Nevertheless, using 
the condition a-y/E'^ <C 1, we can use the approximation 
method to evaluate the quantization integrals based on 
the idea of splitting the whole integration domain [6, a] 
into the intervals [b, r] and [r, a] in each of which only the 
dominating interaction type is taken into account exactly 
while the other integration types are treated as pertur- 
bations. 

We now find a point r that divides the integration do- 
main 6 ^ r ^ a into the domain b ^ r ^ r where the 
Coulomb potential prevails and the domain r ^ r ^ a 
where the long-range potential v{r) prevails. The method 
for choosing such a point is not unique. The most natu- 
ral seems to find a point r where the long-range poten- 
tial v{r) is equal to the Coulomb potential. From this 
requirement, we have r « (ES^/rjiaY^^. 

We can calculate the quantization integrals (for 
a-f/E^ ^ 1) as follows. We calculate integrals (22) by 
expanding the quasimomentum p(r) in a power series in 
the parameters r/a ^ 1 and r/\d\ <C 1 in the domain 
6 < r ^ r and in the small parameters b/r <^ 1 and 
\c\/r <C 1 in the domain r < r < a. Splitting the inte- 
gration interval at the point r « {E^/rjiay^^ therefore 
gives the representation for Ji, 



Ji = aVl - 2A (ji +j2), 



(28) 



where the integrals ji and j2 can be written as follows 
up to terms of the first order in the corresponding small 



parameters r/a, r/\d\ and b/r, \c\/r in the expansions for 
the quasimomentum p{r): 



b 

a 

32 = j s/{a-r){r-d) 



1 — -r + ... 



2ad 



dr, 



1-^ + ... 

2r 



dr. 



(29) 

Calculating the table integrals in (29) and collecting 
terms with like dependence on u, we obtain 



Ji = (jyj-ad{l - 2A) 
— \/ —be arccos 



& + C, ({a-d){c-b) 
log 



2 °V 16ad 

b^ c\ ^ a 4- d 



1 



b-c 

(a - df 



d + a 



X arccos 



d-, 



+ 



(a-hd)(&-Fc) 

(77 



E'^ 



(30) 



We note that when the asymptotic expressions for j\ and 
j2 are added, the result does not contain the parameter 
r. 

To expand the integral J2 in the small parameter 
a"f/&, we represent it as a sum of two terms. 



J2 = - 



k 



2ay/l - 2A 



(31) 



where the integrals j\ and j2 can be written in the forms 



~ ^ 1 } dr 
V-ad J {r+p)^{r - 



32 



^{a - r){r - d) 



b){r-c) 



1 1 
J.2 j.(^j.^qj 



(32) 



dr, 



where p = (/{E + m), q = (E + m)/cr(l - 2A). An 
elementary calculation of the integrals in (32) results in 



arccos 



J2 = - 



k 



E+m 

b-c 



2\a\Vl-2X 



E+m 



\E^ 



(33) 



Adding expansions (30) and (33) and combining terms of 
like orders in a, we obtain the transcendental equation 
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determining the energy spectrum from (9), 



2(7(2A-1) 'V2a-(2A-1) 

log 



- m2 " \^4e(£;2-m2^2 

sgn fc / — rn£ \ / 1 

arccos ( — ^ I = ( + - I tt. 



7 arccos 



where 



?7 = (1 — 2A) arccos (771/772)- 



(34) 



(35) 



Although Eq. (34) is much simpler than "exact" qua- 
siclassical equation (23) for the energy levels, solving it 
still requires numerical calculations. Below, we consider 
scivcral limiting cases where Eq. (34) is simplified and 
can be investigated analytically. 

For the parameter values tr < 0.2 GeV^ and 0.3 < ^ < 
0.8, the condition ^ m is well satisfied for all possible 
values of the level energies -B„^fe of heavy-light mesons. 
If we expand the left-hand side of (34) in m/i? <C 1 
up to third-degree terms, we obtain the transcendental 
equation for -E„^fe: 

[(1 - \)A -X]E'^ + 2mE(l - X)(\A - 1) 



o\k\{\-\r 

iv -r <; log 

Am^ + A[Am^ - 2c7^(l - \)\A = 0, (36) 



-2<7(l-2A)(^7riV + Clog 



where 



arccos 



A = 



(i^) 



V1-2A ' 



iV = n. + l + ?^ + l(7arccos(-^)-a. 



(37) 



Solving this equation by the method of consecutive itera- 
tions, wc obtain the desired expression for the eigenvalues 
En^k in the first approximation (up to terms of the order 
Oiaj/E^)): 

i;WKB(a.) ^^-if^B+[B^ + C [2a{l - 2A) 



a\k\n-X) ,A 
' '1 o + 3^ -h X^A + ttN) 



+ Xfh'^{l-XA) 



1/2 



(38) 



where 



C = (1-A)2A- A- 
B= (1- A)(l- AA)m- 



2aC(l - 2A) 

4ag(2 - 2A) 
E(o) 



and £;(o) = - AV^o- Here, E'-^^ is the zeroth approx- 
imation for the energy on which the quantity En^k de- 
pends rather weakly, and we can set E^^^ w -E'nrfc(OI{=o 
in most cases. 

We have obtained formula (38) for the energy levels 
En^k, which depend nonanalytically on the string tension 
a and which therefore cannot be obtained in the pertur- 
bation theory framework. We mention that for a purely 
scalar confinement (A = 0), formula (38) is simplified to 



j^WKB(as) 



2 

TT 



m + \ m? + a-K log 



a\k\ 



{2E(o)y 



■ttN 



.(39) 



The results of calculating the energy levels E^^^ and 

e"^^^^^"^ based on transcendental equation (23) and 
asymptotic formula (38) together with the exact values of 
Erirk obtained by solving the Dirac equation numerically 
are presented in Table I for rir = 0, 1, 2 and k = ±1, ±2. 
In these calculations, we set the values of ag, A, Vq, mu,d, 
and THs to those used in QCD to describe the states of 
B{bu or bd) and Bs{bs) mesons. As can be seen in Ta- 



ble I, the quasiclassical values E^^^ and E, 



WK13(as) 



Urk 



en- 



sure the respective 1% and 2% accuracies (except the en- 
ergy of states with the radial quantum number = 0, for 
which the accuracy of both formulas is about 8%). The 
accuracy of determining En^k from quasiclassical formula 
(38) is therefore such that the first-order approximation 
usually suffices for practical purposes. 

In order to find the dependence of the energy eigen- 
values En^k on the Coulomb coupling constant ^, we 
solved transcendental equation (23) numerically with the 
following choice of parameters determining the form of 
initial interaction potentials (12): = 0.3, A = 0.3, 
Vo = -0.45 GeV, and rn„,d = 0.33 GeV. The graphs of 
dependences of energy levels on the ratio ^/|fc| are shown 
in Fig. 3, where solid lines indicate the dependence of sev- 
eral lowest levels (n^ = 0) with the given value of k and 
dashed lines correspond to the excited states (n^ = 1). 
As could be expected, as the Coulomb parameter ^ in- 
creases, level energy decreases monotonically and devel- 
ops a square-root singularity as ^ — >■ |A;|. This is a man- 
ifestation of the "falling to center" phenomenon for the 
Dirac equation in composite field (12) with the vector po- 
tential V{r), which has the Coulomb singularity at zero, 
V{r) w Vcoui('') = — ^/r as r — >■ 0. As is known (see 
Sec. Ill), every cut-off of the potential V{r) at small dis- 
tances removes the square-root singularity in the energies 
En^k^ and the curve of the level of En^kiO can then be 
smoothly continued into the domain E < 0. 

It can be seen from Fig. 3 that for the states with the 
same n-r, the energy levels with fc > lie much higher 
than levels with A; < 0. This is the influence of the cen- 
trifugal barrier (for instance, this barrier is absent for 
states with A; = — 1, while it suppresses the probabil- 
ity of the presence of the quark at large distances for 
states with k = +1). These conclusions are completely 
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TABLE I. The results of calculating the level energies E^j^ (based on transcendental equation (23)) and k (based on 

quasiclassical expression (38)) and also the exact values of E„^k calculated at the parameter values as = 0.3, A = 0.3, Vq = — 0.45GeV 
and m„_d = 0.33 GeV, rus = 0.5 GeV (the energies are measured in GeV). 



bu, bd 


bs 


Lj {Ur, k) 


Efi k, 


pWKB 


„WKB(as) 
nrk 




pWKB 


„WKB(as) 

nrk 


Sl/2 


(0, -1) 


0.4327 


0.4408 


0.4729 


0.5248 


0.5322 


0.5623 


(1,-1) 


0.8796 


0.8838 


0.8943 


0.9750 


0.9791 


0.9912 


(2, -1) 


1.1978 


1.2009 


1.2066 


1.2946 


1.2976 


1.3049 


P3/2 


(0, -2) 


0.7355 


0.7373 


0.7504 


0.8376 


0.8392 


0.8460 


(1, -2) 


1.0880 


1.0892 


1.0947 


1.579 


1.590 


1.1927 


(2, -2) 


1.3658 


1.3667 


1.3699 


1.4650 


1.4659 


1.4685 


Pl/2 


(0, 1) 


0.7249 


0.7293 


0.7030 


0.8235 


0.8278 


0.7985 


(1,1) 


1.0701 


1.0733 


1.0594 


1.1696 


1.1728 


1.1572 


(2, 1) 


1.3470 


1.3496 


1.3405 


1.4466 


1.4492 


1.4390 




(0, 2) 


0.9661 


0.9671 


0.9343 


1.0655 


1.0665 


1.0315 


(1,2) 


1.2588 


1.2596 


1.2385 


1.3583 


1.3591 


1.3369 


(2, 2) 


1.5058 


1.5066 


1.4914 


1.6052 


1.6059 


1.5901 



confirmed by numerically solving Dirac system (1) with 

potentials (12); the results of this were presented in [26]. 

We also note that the energies of the lowest levels 
{rir = 0) with fc < reach the zero level {E = 0) at 
the maximum possible value of the Coulomb coupling 
constant ^ = —k (see Figs. 3a-3e). All other states also 
have the singularity of the square-root type at ^ = 
but their energies remain positive. 

The above study of the spectrum of Dirac equation 
in composite field (12) using the WKB approximation 
is of practical interest because calculating integrals in 
quantization condition (9) is much easier in many cases 
than finding exact values of energy levels by numerically 
solving system of radial Dirac equations (1). 



V. THE MASS SPECTRUM OF HEAVY-LIGHT 
QUARK SYSTEMS 

The qualitative picture of forming bound states in a 
Qq system is determined by the presence of the scale pa- 
rameter Aqcd of the confinement of the light antiquark 
Q- ^QCD ^ "^Q) where mg is the mass of the heavy 
quark Q. Under this condition, the heavy quark Q af- 
fects the light quark g as a local static somcc of the color 
(gluon) QCD field. The presence of a small parameter 
AqcD/mg ^ 1 allowed developing powerful means for 
studying QCD in interactions between heavy and light 
quarks. For example, a consistent scheme of the effective 
theory of heavy quarks for hadronic systems with one 
heavy quark {Qq, Qqq) was developed (sec, e.g., [6] and 
the references therein). In the leading term of this the- 
ory (i.e., in the static limit as mg — >■ oo), first, the spin 
of the heavy quark Q splits from the interaction with 
weakly virtual gluons, second, the effective Hamiltonian 
exactly corresponds to the Dirac Hamiltonian of one- 
particle problem (1), and the energy of the spin. orbital 
interaction of the light antiquark q becomes the leading 



term of spin interactions. This is manifested in the ap- 
proximate Isgur-Wise spin symmetry [36] for the heavy 
quark. 

In the leading order in 1/mg, the mass spectrum of 
meson states with one heavy quark is given by the ex- 
pression [6, 37-39] 

MtriQq) = E^^u + ^jEl,-ml+ml, (40) 

where mg and niq are the masses of the heavy quark Q 
and the light quark q constituting the Qq meson. Calcu- 
lating the mass spectrum of Qq mesons therefore reduces 
to consistently calculating the energy eigenvalues of Dirac 
equation (1) in composite field (12) whose source here is 
the heavy quark Q. 

The symmetry properties of Dirac equation (1) drasti- 
cally simplify the problem of classifying states of heavy- 
light mesons. Because the Hamiltonian of Eq. (1) does 
not contain terms describing the interaction of the spin 
of the Q quark with the orbital and spin moments I and 
Sq of the light antiquark, both the spin moment Sq of the 
heavy quark Q and the total moment j = Sq + I oi the 
light antiquark q are two separate integrals of motion. 
This allows classifying the states by the quantum num- 
bers j = ^ , I , . . . of the operator of the total moment of 
the light antiquark g, while the states of the total moment 
of the composite Qq system J = j + Sq are degenerate 
with respect to the orientation of the spin Sq of the heavy 
quark Q. Two almost degenerate states of the compos- 
ite Q^-system with J = j ± 1/2 in the spin symmetry 
approximation [36] therefore correspond to each state of 
the Dirac equation with the given j and with the spatial 
parity P = ( — 1)'+^. Masses of the j^-states of the Qq 
meson are also degenerate with respect to J, and these 
states therefore have identical wave functions. 

The values I = (S states in the quark-antiquark 
model) and j = l/2~ correspond to the ground state 
of the Qq meson. This doublet consists of two states 
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FIG. 3. (a-c) The dependence of the level energies En^k = 
En^k/\Ay on C/I'i^l- Solid lines correspond to the lowest lev- 
els {rir — 0) with the given value of k, and dashed lines 
correspond to the excited states (n^ — 1). The parame- 
ter values are (a) m=0 and A = 0, (b) m = 0.33 GeV and 
A = 0, (c) m = 0.33 GeV, a = 0.18 Gev^ and A = 0.3. Here, 
m*—m/^/a. (d) The dependence of the level energy eo,-i on 
£,/\k\ for different values of the parameter A at m = 0.33 GeV. 



= (0~,1~). In the case I = 1 (the P state in the 
quark model), we have two states with j = 1/2"*" and 
j = 3/2+ and two corresponding doublets = (0+, 1+) 
and = (l+,2+). 

As usual, we introduce a concise notation for the fam- 
ilies of D and Ds mesons: {Dq, D[) are the components 
of the charmed doublet = (0+,l+) with j = 1/2+ 
for nonstrange states (the cu system), {Dig, D'^^) are the 
components of the same doublet for strange states (the 
cs system), and (Di, D2) and {Dgi, DI2) are the com- 
ponents of the doublet = (1+, 2+) with j = 3/2+ for 
the respective nonstrange and strange states. We also use 
the analogous notation system for B and families. 

Above, we did not take the level hyperfine structure 
(HFS) into account, and the proposed potential model 
can predict only the position of the center of masses 
of the HFS multiplet comprising sublevels with different 
moments J = j + Sg. In actual Qq systems, the degener- 
acy of doublet states corresponding to different moments 
J = j ± 1/2 at the given j is broken primarily because of 
the Sq Sq interaction. Therefore, to be able to compare 
our theoretical predictions with experimental data, we 
present the observation values for the centers of masses 
of the HFS multiplets in Tables II-V; these centers of 
masses are calculated by the known formula 



M, 



exp 



E((2J+l)Afj) 
_J 

E(2J+1) 

.7 



(41) 



where Mj is the experimental value of the mass of state 
with the given J. 

Based on these observations, we have tried to describe 
the spectra of masses of low- lying states of the heavy-light 
B{buoTbd), Bs{bs), D{cuor cd), and Ds{cs) mesons con- 
sidering a and A to be universal quantities and setting the 
values of the parameters as and Vq constant in every fam- 
ily of heavy-light mesons allowing them to vary slightly 
only when passing from one family to another. All the 
parameters a, A, as, and Vq of potential model (12) were 
determined by fitting the known data for the mass spec- 
tra of pseudoscalar D and B mesons. The found values 
of the parameters are consequently used below in other 
applications in the framework of our approach, for exam- 
ple, when describing the spectra of the strange Dg and 
Bs mesons. 

We use only one a priori restriction: the value of the 
coefhcient A of mixing between the vector and scalar long- 
range potentials Vconiir) and Scon{{T) must lie in the in- 
terval ^ A < 1/2 for the EP U{r,E) of interaction 
model (12) to be a confining-type potential. The value 
of the parameter A was obtained by fitting experimental 
data [40, 41] on the fine structure of P-wave levels in D 
and B mesons. It was established that the fine struc- 
ture of the P-wave states in the heavy-light {D, Dg, B, 
and Bg) mesons is primarily sensitive to the choice of the 
mixing coefficient A and to the value of the strong cou- 
pling constant as . Comparing the results of calculations 
based on formulas (23) and (40) with the experimental 
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data [40, 41], wc find that the best agreement is reached 
at A = 0.3 and for the parameter choices 

a = 0.18GeV^ ^^(cMorcfi) = 0.386, a,{buoxhd) = 0.3, 
Vo{cu or cJ) = -375 MeV, Vo{hu or hi) = -450 MeV. 

For the masses of w, d, s, c, and b quarks, we used their 
constituent masses rriu.d = 330 MeV, = 500 McV, 
mc = 1550 MeV, and nib = 4880 MeV. When calculat- 
ing the mass spectrum, we neglected electromagnetic in- 
teraction and the difference of the masses of u and d 
quarks, therefore considering the particles Z)+, , D'^ , 
and -D", for example, to be the same state of the Qq sys- 
tem, = 0~. Correspondingly, we do not distinguish 
between the interaction parameters cr. A, ofg, and Vo for 
these particles. The mass spectra of D and Dg mesons 
calculated in this approximation and by means of numer- 
ical solutions of system (1) are presented in Tables II and 
III. 

TABLE II. The mass spectrum and the mean radii of D mesons 

obtained in the WKB approximation and numerically for potentials 
(12) (masses are expressed in MeV and the mean radii are expressed 
in Fm). 



Lj (rir, k) 


Mnum 




Mexp 


(r) num 


(r)(43) 


Sl/2 


(0, -1) 


1989.1 


2001.5 


1971.1 


0.472 


0.402 


(1,-1) 


2624.5 


2632.3 


<2637 


0.684 


0.664 


P3/2 


(0, -2) 


2440.1 


2443.2 


2447.3 


0.678 


0.632 


(1, -2) 


2979.7 


2981.9 




0.856 


0.833 


Pl/2 


(0, 1) 


2395.2 


2403.7 


2407.8 


0.513 


0.568 


(1, 1) 


2926.8 


2933.4 




0.770 


0.788 



The agreement between the numerical result (result of 
WKB approximation) and experiment is less than 1% 
(1.5%), except for the masses of states P3/2 a-nd P1/2 of 
the cs system for which the mismatch depends on the in- 
terpretation of the r>si(2536)^ meson and is 3.1% (3.5%) 
if we consider it to be the vector state = 1+ belonging 
to the doublet j = 3/2+ or 0.4% (0.3%) if we consider 
it to be the state = 1+ of the doublet j = 1/2+. 
There is a rather broad spectrum of opinions concerning 
the identification of the states P3/2 and P1/2 of the me- 
son with the quark content cs (see, e.g., [42-52]). For 
example, the state = 2+ of a relatively narrow dou- 
blet j = 3/2+ is related to Z?s2(2573), while the vector 
state ~ 1+ belonging to the doublet j = 3/2+ is 
generally related to the isotopic singlet £)si(2536)+ me- 
son with the mass 2535.35±0.34±0.5 MeV (values (I) for 
Mcxp in Table III) [42-50]. On the other hand, the state 
£'si(2536)+ was associated with the state = 1+ of 
the wide doublet j = 1/2+ in [51, 52] (values (II) in Ta- 
ble III). Therefore, a reliable experimental identification 
of this state is still lacking. We note that our calculations 
agree better with the second possibility. 

For bu and bs systems, we obtained a good agree- 
ment of our results with the experimental data for the 
ground state with j = l/2~ and for the P state with 
j = 3/2+ (see Tables IV and V). For states in the dou- 
blet j = 1/2+ , we have only theoretical predictions of 



other authors. For the bu system, our restilts agree with 
the data obtained in [53], and a remarkable agreement 
with the results in [47, 48, 54] was obtained for the bs 
system. 

In the leading approximation (in l/mg), the wave 
functions and excitation energies of the strange quark 
in the field of a heavy c or 6 quark reproduce the corre- 
sponding characteristics of heavy-light mesons with light 
u and d quarks with high accuracy. Therefore, up to 
an additive upward shift of masses on the value of the 
current mass of the strange quark 

m,, w M[Ds] - M[D] « M[Bs] - M[B] « 0.1 GeV 

the level systems for Dg and Bs mesons coincides with 
the respective level systems for D and B mesons if we 
do not take the level splitting depending on the spin of 
the heavy quark into account. Further, the spin-orbital 
splitting of lower states of Dg and Bg mesons for the 
levels P3/2 and P1/2 is 35% larger than that of the D 
and B mesons. 

Not only the spectrum of bound systems, all other ob- 
servable characteristics of heavy-light mesons can be cal- 
culated in the framework of the quasiclassical approach 
under consideration. For example, an important meson 
characteristics is the mean radius (r), which determines 
the radius of the light quark orbit in a definite state \nrk) 
in the case of hydrogen-like quark systems. Wc first ob- 
tain general formulas expressing the means of type (r™) 
(i.e., the moments of the probability distribution density) 
in terms of quasiclassical asymptotic expressions for so- 
lutions of the Dirac equation. Using the standard proce- 
dure, we obtain the known quasiclassical formula 

00 00 
{rn= jx^r"'xdr = j {\F {r)\^ + \G{rf)r^dr 



2 [E-V{r) 
p{r) 



r'^dr, 



(42) 



where the period T of radial oscillations of the classical 
relativistic particle is given by the formula T = 2 J^' {E- 
V{r))/p{r)dr[lQ]. 

All the integrals in (42) can be expressed in terms of 
complete elliptic integrals (A.l). In particular, the mean 
radius of the bound state is 



{r) = 



4 [niF (x) + n2E (x) + nsH {v, x)\ 
T^{a-c){b-d){l-2\)\a\ ' 

^ 4 [mF (x) + n^E (x) + noR {v, x)\ 
^{a-c){b-d){l~2X)\a\ 



(43) 



(44) 



where the quantities n, (i = 1, . . . , 6) are defined in Ap- 
pendix. The calculation results for (r), according to for- 
mulas (43) and (44) for different states of D, Ds, B, 
and Bs mesons are presented in the last columns in Ta- 
bles II-V. We see that the quasiclassical approximation 
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TABLE III. The mass spectrum and the mean radii of Ds mesons obtained in the WKB approximation and numerically for potentials 
(12) (masses are expressed in MeV and the mean radii are expressed in Fm). 





Mnum 


MwKB 


Mexp 


(''")num 


{r)(43) 


Sl/2 


(0. -1) 


2057.2 


2069.0 


2072 


0.416 


0.359 


(1. -i) 


2729.1 


27:i7.1 




0.616 


().()28 


P3/2 


(0, -2) 


2550.1 


2552.1 


2559.2 (I) 


2530.7 (II) 


0.625 


0.588 


(1, -2) 


3105.2 


3107.2 






0.814 


0.795 


Pl/2 


(0, 1) 


2499.7 


2508.5 


2423.8 (I) 


2480.9 (II) 


0.504 


0.536 


(1, 1) 


3051.7 


3058.5 






0.739 


0.756 



TABLE IV. The mass spectrum and the mean radii of B mesons 

obtained in the WKB approximation and numerically for potentials 
(12) (masses are expressed in MeV and the mean radii are expressed 
in Fm). 



LA 


rir, k) 


A^num 


A^WKB 


Mexp 




(r)(43) 


Sl/2 


(0, -1) 


5320.7 


5329.5 


5313.5 


0.516 


0.448 




(1,-1) 


5827.3 


5832.2 




0.728 


0.708 


P3/2 


(0, -2) 


5659.6 


5661.6 


<5698 


0.711 


0.666 




(1, -2) 


6076.9 


6078.4 




0.888 


0.865 




(0, 1) 


5647.4 


5652.4 


5751.6 [47] 


0.577 


0.612 


Pl/2 








5624 [53] 








(1, 1) 


6055.1 


6059.0 




0.812 


0.829 



TABLE V. The mass spectrum and the mean radii of Bs mesons 
obtained in the WKB approximation and numerically for potentials 
(12) (masses are expressed in MeV and the mean radii are expressed 
in Fm). 





Tlr, k) 


Afnum 


*^WKB 


Mexp 


('")num 


(r)(43) 


Sl/2 


(0, -1) 


5407.4 


5415.6 


5404.8 


0.457 


0.404 




(1,-1) 


5926.3 


5931.2 




0.688 


0.671 


P3/2 


(0, -2) 


5763.7 


5765.6 


<5853 


0.656 


0.619 




(1, -2) 


6185.5 


6186.8 




0.845 


0.826 










5751.8 [47] 














5753.3 [48] 






Pl/2 


(0, 1) 


5747.2 


5752.2 


5700.5 [53] 


0.547 


0.575 








5755.0 [54] 
5790.3 [55] 








(1, 1) 


6162.8 


6166.8 




0.779 


0.795 



well describes the numerical simulation results (r)num 
and ensures an accuracy up to 3% (except the ground 
state). Calculations demonstrate that the mean radius 
of the Qq system increases monotonically as the energy 
increases. 

In addition to the "exact" quasiclassical formulas (43) 
and (44), it is desirable to find approximate analytic 
expressions for the quantities (r) and T. We already 
addressed an analogous problem in the preceding sec- 
tion when constructing asymptotic approximations for 
the quantization integrals. 

If the condition cr/^m^ <C 1 is satisfied in the spectral 
domain E < m, then the only essential contribution to 



the integral determining the mean radius (r) comes from 
the domain of the integration variable r where the long- 
range potential v{r) can be considered a small perturba- 
tion. Neglecting this potential, we obtain the expressions 
for the mean radius and the period in the zeroth approx- 
imation: 

A more accurate expression (than (45)) for the mean ra- 
dius can be obtained if we use the exact solutions of 
Dirac system (I) in the Coulomb field in the integral 

/(r(l^(^)l^ + \Gir)f)rdr = (r) [20]. The resulting ex- 
pression for the mean radius of the hydrogen-like system 
becomes 

and it coincides with (45) at large values of the radial 

quantum number rir- 

This simple approximation ensures an amazingly good 
accuracy for deeply lying levels (but, of course, not at 
E = 0). For example, for the first three terms IS1/2, 
lPi/2, and 2S1/2 of the b quark {mi, = 4.88 GeV), we 
obtain the respective values (r) = 0.153 Fm, 0.501 Fm, 
and 0.609 Fm from (46), while the exact calculation (the 
numerical solution of the Dirac equation with potentials 
(12) at ^ = 0.4, A = 0.3, Vo = -0.45 GeV, and a = 
0.18 GeV^) yields the respective values (r) = 0.153 Fm, 
0.493 Fm, and 0.600 Fm. Our approximation therefore 
ensures a high accuracy in the case of heavy quarks. 

Unfortunately, the domain of applicability of such an 
approximation is restricted by the condition a/S,m^ ^ 1. 
Because the problem of a size of a bound state of the 
Qq system is important, we consider it from the quan- 
titative standpoint. We use the fact that the condition 
aj/E"^ ^ 1 is satisfic;cl for all typical values of the pa- 
rameters ^ and a of heavy-light quarks in the spectrum 
domain E > fh under investigation. In this case, the light 
quark motion is mainly determined by the linear poten- 
tial, and the Coulomb interaction can be considered a 
perturbation. In some cases, the zeroth approximation 
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suffices for calculating (r) and T, 

2 ( ^ 3Xr]i 



(r) 



Tcr2(l-2A) 

irccos ?^ r ~ 

Em + 



V1-2A 



2(1 - 2A) 



2 V 1 - 2A 



(7(1 - 2A) 



-A 



+ (1-A)ryry2 



,(47) 



(48) 



where the quantities rji and 772 are determined in (25), 
and the quantity 77 is determined in (35). For example, 
for the first three terms IS1/2, lPi/2j and 2Si/2 of the B 
meson {mi, = 4.88 GeV and m„ = 0.33 GeV), we obtain 
the respective quantities (r) = 0.381 Fm, 0.576 Fm, and 
0.681 Fm in approximation (47), (48), and the calcula- 
tion using "exact" formulas (43) and (44) (at ag = 0.3, 
A = 0.3, Vo = -0.45 GeV, a = 0.18 Gev^) yields the re- 
spective values 0.448 Fm, 0.612 Fm, and 0.708 Fm. This 
approximation therefore ensures an acceptable accuracy 
for calculating the mean radii of the Qq mesons. 

In contrast to the masses of bound states of the heavy- 
light mesons, the wave functions of the mixed mesons 
at the zero, which determine the leptonic constants and 
normalizations of cross-sections of formation D- and B- 
mesons, are more sensitive to global properties of initial 
potentials S{r) and V{r). We now turn to deriving these 
wave functions. 



VI. ASYMPTOTIC COEFFICIENTS OF A 
WAVE FUNCTION 



Asymptotic coefficients Cf.g a-t the zero and Ap^c at 
the infinity are the characteristic parameters of a wave 
function. We shall give simple analytical approximations 
for these coefficients which describe the results of numer- 
ical calculations quite well. 

First of all we consider the construction rules of asymp- 
totic expansions of solutions of the Dirac equation at zero 
(r — ^ 0) along with more standard expansions of solu- 
tions, when r — >■ 00. For the considered potentials (12) 
we have 

F,G = CF,Gr'' + ...,r^O, Cf/Cg = {k-j)/C, (49) 

and for wave functions of the discrete spectrum (0^ A < 

00 

1/2) the normalization condition J (F^ + G^) rfr 1 is 



satisfied. Values of ^^(0), 0^(0) (or, more precisely, 
of Cpq) define the probability to discover particles at 
small distances from one another and are of consider- 
able physical interest especially in the case of systems, in 
which interactions of two various types (for example, the 
Coulomb interaction and long-range one v(r)) exist. 

In the classically forbidden range < r < ro the wave 
function of oscillating type is changed by the solution 
decreasing exponentially with increasing r (see Fig. 1). 
Matching the WKB-solutions of the Dirac equation on 



both sides of the turning point tq, wc obtain the quasi- 
classical expressions for the radial wave functions F{r) 
and G{r) in the classically forbidden region < r < tq: 



F{r) = i-iy 



E-V + m + S' 

q{r) 

ra 



1/2 



X exp 



q 1 dr 

qr 



,(50) 



G(r) =sgn A; (-1)"-^ 



E-V-m-S 



q{r) 



1/2 



X exp 



q 1 dr 

qr 



.(51) 



All integrals in (50) and (51) are expressed through the 
quite complicated combination of the elliptic integrals. 
But in the cases Er <fh and Er>m they can be calcu- 
lated through elementary functions, using the relations 
ajijv? -C 1 and a^/E'^ <C 1 to expand the integrands 
into power series. 

Let us first investigate the asymptotic behavior of 

the quasiclassical solutions (50), (51) at r — >■ for the 
lower levels {E < rh) which are basically defined by the 
Coulomb potential [al^rv? <^ 1). Note that the larger 
the Coulomb parameter ^, the smaller is the essential 
potential v{r) at small distances. Before the evaluation 
of the asymptotic coefficients Ci?,G by means of formu- 
las (50), (51) it is necessary to expand the quasiclassical 
momentum p{r) in potential v{r). Then, using the tech- 
nique of evaluation of the phase integrals from Sec. IV 
and proceeding in (50), (51) to the limit r — ^ 0, we obtain 
in zeroth approximation the expressions for the asymp- 
totic coefficients at zero: 



Oo{\k\-i) 




^{jm+\k\Eo) 
fc — 7 



Cf 
Cg 



(52) 



Here 6q = \J {Eq ky — {m , and the period of ra- 
dial oscillations T is given by the previous formula (45). 
When wc derive the expression (52) we use the qua- 
siclassical requirement of the normalization (5). Solv- 
ing the Dirac equation (1) at small distances (in range 
< r < c) one can neglect the term with the linear po- 
tential (cr = 0). Having used the asymptotic behavior 
of the normalized radial functions F{r) and G{r) of the 
relativistic Coulomb problem [20] at r — )■ and the rela- 
tions (49), we find the more exact (than (52)) expression 
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for Cf 



Up — 



7+1/2 



r(27 + i) 



Ho 



Mo 



k 



1/2 



.(53) 



Here n[. = n,. + (1 + sgnfc)/2. The formulas (52) and 
(53) differ one from another within an error between the 
Stirhng formula 



n! = V27rexp{(n + 1/2) logn - n} [1 + 0{n )], n ^ oo 

and the F-function. ^ 

For states with Er > m, when the requirement 
a-y/E^ ^ 1 is satisfied, the Coulomb potential is essential 
only in the range of small distances and can be consid- 
ered as a small perturbation in the basic range of particle 
localization (i.e. in classically allowed range c < r < b). 
This gives the possibility to exclude (cr = 0) the linear 
part of the potential v{r) from the quasiclassical momen- 
tum p(r) when evaluating the integrals in exponents (50), 
(51). Then the asymptotic behavior (at r — > 0) of radial 
wave functions F{r) and G{r) obtained in this way allows 
to determine the asymptotic coefficients: 



\Cf\ 




H\k\-ll 
C(7m + \k\Er 



X exp 



^Er 



: arccos 



^Er 



Ei 



, (54) 



where quantity is defined in (24), and energy is 
given by the formula (40) . Characteristic feature of the 
considered case is the fact that in the integral (5), which 
defines a period of radial oscillations T, only the range of 
values of the integration variable r, where the Coulomb 
potential can be considered a perturbation, is essential. 
By neglecting the Coulomb interaction, we arrive at the 
previous expression (48). 

The asymptotic coefficients Ap, Aq of radial wave 
functions at infinity are important physical parameters 
of bound states as well. Along with the coefficients Cp, 
Co at zero (49), the asymptotic coefficients Ap^c are con- 
tinually encountered in quantum mechanics [22], atomic 
and nuclear physics [56, 57], in the converse problem of 
quantum scattering theory [58, 59] etc. For the poten- 
tials (12) the quantities Af,g are related to asymptotic 
behaviors of the normalized radial wave functions by the 
relations 



F,G = Ap^cr^exp 



V1-2A 



where ar 

2 
V2 



oo 



a > 
AC 



j55) 
7 = 



0, Ap = -Vi^^Ag, 
, and the parameter A has val- 



2(1 - 2A)3/2a VI - 2A 
ues in the range ^ A < 1/2 

In the below-barrier range r > ri = b far from the turn- 
ing point ri = b under the requirements aj^TV? <SL 1 and 
E < fh, after the evaluation of integrals the quasiclassical 
solutions (6) are of the form of decreasing exponents 




Vl - 2A(?o + (1 - 2A)crr + r/i + ^^(1 - 2A)(t//^2 



C- 1 \ 1/4 / t? \ 



yr^2A[M + ^{\^J? + 2r]iE)/ii3] + ryi + ^E(l - 2\)a/iJ? 



exp 



2(l-2A)a 



2/^2 



(56) 



where qq — (1 — 2A)(r — c)(r — d). The estimates 
show that by the satisfaction of the requirements 
uj^rr? <C 1 and E < m there is quite long range of dis- 
tances r which are much larger than size of the Coulomb 
hydrogen-like system (r » (r), see (45) or (46) and much 
smaller than the distance r sa {E^/rjiay/"^ at which the 
Coulomb interaction becomes quantitatively comparable 
with the long-range interaction. In this range as the wave 
functions of zeroth approximation it is natural to take the 
unperturbed radial functions F and G of the relativistic 
Coulomb problem, and the potential v{r) can be consid- 
ered as a small perturbation. Neglecting it, we arrive at 



the following quasiclassical expressions for F and G 



m + Ea \ WKB(as) ^ -anr 




m + Eq 



— fc/io 
kno 



1/4 



£,Eo - 7/xo 



2^gr^ ^0 
~0o 



(57) 



Equating (57) to the asymptotic (at r — > oo) represen- 
tation of solutions of the Dirac equation in the Coulomb 
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field [20] 




Acr t'o e 



(58) 



we obtain the explicit expression for the period of radial 
oscillations of the classical relativistic particle: 



T = 



1 



1/2 



/ 2ei4 



V ^0 



(59) 



Here Ac is the asymptotic coefficient of the Dirac radial 
wave functions in the Coulomb potential: 



\Ar 



-,1/2 



2^m2r(27 + < + l)<! 



(2mo) ^0 ■ 



(60) 



Comparison of (57) and (58) shows that their exponential 
and power factors are the same, however, the asymptotic 

coefficients differ within an error be- 

tween the Stirling formula and F-function. 

The formula (57) is obtained by neglecting the linear 
part of the potential v{r). This approximation was ar- 
gued above by nicians of the circumstance that under the 
quasiclassical requirements {a/S,fh^ "^1) there is a range 
of distances (r) <C r <C r, in which the distortion of a 
wave function, caused by action of the linear part of the 
potential v{r), can still be neglected and there is the law 
of decreasing radial wave functions (57) that is charac- 
teristic for the relativistic Coulomb problem. Change of 
the law of decreasement (57) of functions F{r) and G{r) 
to (55) at r » f appears because in EP U{r, E) we have 
taken into account the quadratic (in or) terms which in- 
crease with increasing r more rapidly than others and so 
play a role of the perturbation which destroys the asymp- 
totic regime (57). As a result of such an account, using 
the quasiclassical approximation (56) for the normalized 
radial wave functions F and G at large r, we obtain the 
following expression for the asymptotic coefficient at in- 
finity 



2mo 



^c(l-2A)^+V^(^^^^° + ^^")"''° 

(yr^2A/x-r?i)2 



Xcr ' 



.7 / /^O 



MO 



exp 



4(1 - 2A)3/V 



(61) 



We now proceed to the other limiting case a^jE"^ <^ 1 
when the centrifugal potential 7^/2mr^ does not play an 
essential role at large distances and can be omitted in the 
quasiclassical momentum p(r) = iq{r). Having expanded 
the quantity q{r) = \p{r)\ in powers of the Coulomb po- 
tential and calculated the integrals in exponents in (6), 
under the requirement a^jE^ <C 1, in the asymptotic 



domain r — > oo we arrive at formulas of type of (55) for 
F and G, in which 



T 



4(1 - 2A)3/V 



, ~ , arccos 

2v'l - 2\{E^ - m2) ^£;2_^2 V ^ 



(62) 



and the period T is determined by the previous expres- 
sion (48). 



VII. ENERGY SPECTRUM OF THE MASSLESS 
FERMION IN THE EXTERNAL SCALAR FIELD 
OF THE FUNNEL TYPE 

In this section, we shall study energy spectrum of the 
Dirac equation (1) for massless fermion (m = 0) in the 
external scalar field of the form 

S{r) = -^+ar, a> 0; V{r) = 0. (63) 

Physical details of the considered model (63) can be 
found in [16]. 

For the potential (63) and particle with zero mass 
the quasiclassical quantization condition (9) becomes the 
transcendental equation 



E^ + 2a(f - 7) 



+ 



4cr 

n(t/l,x*) 



a{b + ajn 
1 



2b 



+ 



b + P+ b + P. 



b^-Pl 



= n, 



1 

+ 2- 
(64) 



Here K(x*) and Il(i^*,x*) are the complete elliptic inte- 
grals of the first and third kind, and the following nota- 
tions are used: 



lE^ + 2a{i' - 7) 



\/£2 + 2c7(e + 7)' ± ^ P±-ro' 

P± = ^{-E±./WTW^). 



where the turning points ro = b and ri = a are deter- 
mined by the equations 

b,a= ^^fE^+2^^^^^/WTW^W^f^, 



7 = a/A;^ + ^'2- In the case a — > the equation (64) 
can be solved in an explicit form. At small values of 
parameter a (namely at cr < 0.2 GeV^) the condition 
E^ aj for all possible values of level energy En^k is 
satisfied very well. In this case the formula given above 
become appreciably simpler and the equation (64) for 
a quasiclassical spectrum ultimately assumes the rather 
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TABLE VI. The eigenvalues E„^h of the massless Dirac equation 
with the scalar potential (63) for ^ = 0.4, cr = 0.18 GeV^ 



States 


-B™™, GcV 




E^lk, GcV 


En%, GcV 


rir 


k 


(1) 


(64) 


(65) 


(66) 





-1 


0.5568 


0.5581 


0.6530 


0.6568 


1 


-1 


1.0293 


1.0298 


1.0489 


1.0489 


2 


-1 


1.3379 


1.3382 


1.3471 


1.3470 


3 


-1 


1.5862 


1.5863 


1.5917 


1.5917 


4 


-1 


1.7999 


1.8000 


1.8037 


1.8037 





-2 


0.8217 


0.8220 


0.8561 


0.8554 


1 


-2 


1.1898 


1.1899 


1.2003 


1.2001 


2 


-2 


1.4649 


1.4650 


1.4704 


1.4703 


3 


-2 


1.6949 


1.6949 


1.6983 


1.6983 


4 


-2 


1.8966 


1.8966 


1.8990 


1.8990 





1 


0.9340 


0.9335 


0.9054 


0.9055 


1 


1 


1.2568 


1.2566 


1.2457 


1.2457 


2 


1 


1.5143 


1.5141 


1.5079 


1.5097 


3 


1 


1.7346 


1.7345 


1.7304 


1.7304 


4 


1 


1.9302 


1.9302 


1.9272 


1.9272 





2 


1.0996 


1.0994 


1.0857 


1.0853 


1 


2 


1.3846 


1.3845 


1.3781 


1.3780 


2 


2 


1.6217 


1.6216 


1.6177 


1.6177 


3 


2 


1.8289 


1.8288 


1.8261 


1.8261 


4 


2 


2.0152 


2.0152 


2.0132 


2.0132 



simple form 



E 



2(7 



ak 



- log - 1 - i? 



+0 



(77 



(65) 



where cr > anti 

iV' = 2nr + 1 + 7 - ^' + 2B sgn fc, 



B = — arctan . 

TT 



7-C 



7'^ 



/^e: 

fc2 |fc| 



In that case, when the Coulomb-like term in the poten- 
tial (63) is absent (^' = 0), the equation (65) coincides 
exactly with the quasiclassical quantization rule for en- 
ergy levels in the scalar well U (r, E) generated by linear 
confining interaction [15] S{r) = ar, V{r) = 0. 

The equation (65) for En^k is easy for solving numer- 
ically. Comparison of results of such calculations En^h 
with exact values [26] (see tab. VI), obtained by numer- 
ical integration of the Dirac system (1), shows that the 
quasiclassical equation (65) provides acceptable accuracy 
of calculation of the energy spectrum: even for the lower 
states with ^ 1 the error (65) in determination of 
Eurk does not exceed 5% and rapidly decreases with in- 
creasing rir- Along with the direct numerical solution of 
transcendental equations (64) and (65) it is worthwhile to 
construct (by means of some simplifications or approxi- 
mations) the approximate analytical expressions for level 
energy that would allow without difficulty to trace depen- 
dence of En^k on quantum numbers rir, k and parameters 



of interaction model (63). Solving (65) by the method of 
iterations we arrive at the analytical expression for en- 
ergy 

s^^, = ^^±{N' +[{N' -a'' 

-|l/2~ 



/(7 



„, , 1 A 8(iv'-e') 

+2A: ( - hog 



R 



■ (66) 



The positive sign of the root corresponds to energy of 
a particle, and the negative one corresponds to the an- 
tiparticle energy taken with a minus sign. 



VIII. SUMMARY 

The most important results of investigation performed 
can be summarized as follows: 

1. The relativistic potential quark model of Qq- 
mesons in which the light quark motion is described 
by the Dirac equation with a scalar-vector inter- 
action and the heavy quark is considered a local 
source of the gluon field is constructed. The effec- 
tive interquark interaction is described by a com- 
bination of the perturbative one-gluon exchange 
potential Vcon\{T) = —^/r and the long-range 
Lorentz-scalar and Lorentz-vector linear potentials 
SLr.{r) = il-X){ar + Vo) and Vi.,.(r) = X{ar + Vo). 
It is established that the quark confinement al- 
ways arises when the Lorentz-scalar part 5i.r. of 
the long-range interquark interaction prevails over 
the Lorentz-vector one VI. r.. 

2. The new asymptotic method of calculation of quan- 
tization integrals that based on splitting an inte- 
gration interval into two segments in each of which 
only the dominating interaction type is taken into 
account exactly while the other interaction types 
are treated as perturbations is elaborated. Approx- 
imative analytical expressions for energy spectrum 
of heavy-light mesons obtained within quasiclassi- 
cal approach at u^jE? <C 1 are asymptotically ex- 
act in the limit n,. — > 00 and ensure a high accuracy 
of calculations even for states with the radial quan- 
tum number ^ 1. 

3. In the framework of the considered model we have 

obtained the satisfactory description of the mass 
spectrum of £)-, Dg-, B-, and iJ^-mesons. We show 
that the fine structure of P-wave states in heavy- 
light mesons is primarily sensitive to the choice of 
two parameters: the strong-coupling constant 
and the coefficient A of mixing of the long-range 
scalar and vector potentials S'i.r.(r) and 14.r.(?')- 
The best agreement between the theoretical predic- 
tions and experimental data exists when the mixing 
coefficient A = 0.3. 

4. Using WKB method the convenient analytical for- 
mulas for asymptotic coefficients of wave function 
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at zero and infinity and mean radii of the Qq 
mesons. 



r dr 



R Via-c) {b-d) 



APPENDIX 

We consider the quantization integral Ji. We rewrite 
the expression for Ji in (22) as the sum of integrals 

Ji = + /i9o + 5^1 + /92 + 93), 



R 



dr, n = -1,0, 1,2,3, .... 



where the quantities /, g, h, and I arc determined in (20) 
and the quantity i?(r) is determined in (21). After the 
standard change of the integration variable [33] 

b{a — c) — c(a — b) sin^ if 
a — c — (a — b) sin^ 

the integrals Sj„ are expressed in terms of the complete 
elliptic integrals of the first, second, and third kind, which 
are written in the conventional notation [34] as 



7v/2 



Fix) 



d(p 

7r/2 



E 



tt/2 

(X)= J A dip, 



(A.1) 



n(.,x) = /(^ 



df 



(1 — i^sin 99) A 



+c{b-c)U{iy,x)+{b-cfT2(l,iy,x) , (A.5) 



«3 



•dr 



R V(«-c) {b-d) 



c'F{x) 



+3c2 (6 - c) n (ly, x) + 3c (6 - cf T2 (f , x) 



+ (b-cfn{l,y,x) 



The integrals of the form 



Tn i^P, V, x) 



dip 



(1 — vsiv? ip) A 



are calculated using the recurrence relation 



Tn-3 — 



'V^ sin Lp cos ipl\ 



(2n-5)x2 [(l-,.sin2<^) 



(A.6) 



A = \Jl-X^ sin^ ip, V = 



a — b 



a — c 



ic-d) 



{b-d)- 



We thus obtain the representations for ...,33: 



+2(n - 2) [3x' - v{\ + x')] T^-i 



-(2n - 3) [x' (3 -2v)^v{v- 2)] T„_i 



3-1 



dr 



rR ^{a - c) (6 - d) be 



bF{x) 



+2{n-l){x^-y){l-u)T,, 



9o 



-(6-c)n(^:.,x) 

a 

-I 



R y (a - c) (6 - d) 



(A.2) 



^(X), (A.3) 



Wc analogously find the integrals in the expression for 
J2 in (22): 



dr 



I {r- 



X±)R y(a-c) (6-d)(6-A±)(A±-c) 



3i 



a 

I 



■dr 



R ^(a-c) {b-d) 



[cF{x) 



{X±-b)F{x)-{b-c)U{^^-f^y,X 



, (A.7) 



+ {b-c)U{u,x)], 



After expressions (A.2)-(A.7) are substituted in inte- 
(A.4) grals (22), quantization condition (9) becomes transcen- 
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dental equation (23), where 



quantities (i = 1, . . . , 6) in formulas (43) and (44): 



^± = ^^1^, 5R = (1 - i.) (x^ - z.) , 



8(1-:/) 



+ 2 + c^f + cg + h + l/c) 



(b-c) 

K = (3 - 2u) + v{v - 2) 

V 



No 



3(6-c)H 



3 (6 - c) H2 23? , 2 „ . N 



+ (6-c)((l + x')!^-3x')+K(/ + 3c) 

^^--(6^^' iV. = [(^-A,)(A,-c): 
7V6 = [(6-A_)(A_-c)]-i, 

2 i^ + TO \ 

N7 = 



(A+-c)(A_-c) I 2(l-2A)cr 



rii = E I c 



,2 ~ "-y 1 Xa i - ~ 



2(1-^) 



+ 



{b-cf 



X - 



3K 



v{b- cf 

m 



E-2,Xa 



(6 - c) H 
4JR 



na = (6 - c) 
- Act 3c' 



E[2c + 



2 (b-C)^ 

25i 



(& - c) H 
2SR 

3c- 



(b-c)(3x'-i^(l + x')) ^ 3(6-c)H 



714 = c£ - Ao- c^ - ^ 



2(1-1/) 



4JR 



Xau (b — cY 

m 



We analogously find the integrals appearing when cal- 
culating the mean radii by formula (42). We present the 



ne = {b- c) 



E - Act 2c + 



(6 - c) H 
23? 
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